Optimal inequalities for state-independent contextuality 
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Contextuality is a natural generalization of nonlocality which does not need composite systems 
or spacelike separation and offers a wider spectrum of interesting phenomena. Most notably, in 
quantum mechanics there exist scenarios where the contextual behavior is independent of the quan- 
tum state. We show that the quest for an optimal inequality separating quantum from classical 
noncontextual correlations in an state-independent manner admits an exact solution, as it can be 
formulated as a linear program. We introduce the noncontextuality polytope as a generalization of 
the locality polytope, and apply our method to identify two different tight optimal inequalities for 
the most fundamental quantum scenario with state- independent contextuality. 
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Introduction. — The investigation of the operational dif- 
ferences between quantum mechanics and classical me- 
chanics resulted 1964 in the discovery of Bell's inequali- 
ties Such inequalities constrain the correlations ob- 
tained from spacelike-separated measurements and are 
satisfied by any local hidden variable (HV) model but are 
violated by quantum mechanics. For every measurement 
scenario, there exists a minimal set of inequalities, called 
tight Bell inequalities, which provide also sufficient con- 
ditions: If all tight inequalities are satisfied, then there 
exists a local HV model reproducing the corresponding 
set of correlations [2|, |3[ . 

Mathematically speaking, each tight Bell inequality 
corresponds to a facet of the locality polytope [3|. This 
means that it is an (p — l)-dimensional face of the p- 
dimensional polytope obtained as a convex hull of the 
vectors representing local assignments to the results of 
the considered measurements. Such a polytope gives all 
classical probabilities associated to a local model for a 
given measurement scenario, and its facets give precisely 
the boundaries of the polytope. In this sense, tight in- 
equalities separate classical from nonclassical correlations 
perfectly. 

Similarly, noncontextuality inequalities 4-6| are con- 
straints on the correlations among the results of compati- 
ble observables, which are satisfied by any noncontextual 
HV model. While the violation of Bell inequalities reveals 
nonlocality, the violation of noncontextuality inequalities 
reveals contextuality 0, H[ , which is a natural generaliza- 
tion of nonlocality privileging neither composite systems 
(among other physical systems), nor spacelike-separated 
measurements (among other compatible measurements), 
nor entangled states (among other quantum states). 

All Bell inequalities are noncontextuality inequalities, 
but there are two features of noncontextuality inequal- 
ities which are absent in Bell inequalities. One is that 
noncontextuality inequalities may be violated by simple 
quantum systems such as single qutrits Q. These vio- 
lations have recently been experimentally observed with 



photons [jj. The other is that the violation can be inde- 
pendent of the quantum state of the systems thus it 
reveals state-independent contextuality (SIC). The latter 
has been demonstrated recently with ququar ts ( four-level 
quantum systems) using ions (lCf . photons and nu- 
clear magnetic resonance [l2j . 

The notion of tightness naturally applies also to non- 
contextuality inequalities. Tight noncontextuality in- 
equalities are the facets of the correlation polytope of 
compatible observables as we will explain below. Com- 
pared with the locality polytope, the difference is in the 
notion of compatibility, since now one no longer consid- 
ers only collections of spacelike-separated measurements, 
but admits more generally the measurement of a context, 
i.e., a collection of mutually compatible measurements. 
For a given contextuality scenario, the corresponding set 
of tight inequalities gives necessary and sufficient condi- 
tions for the existence of a noncontextual model. 

For example, the three inequalities with state- 
independent violation introduced in Ref. are all tight. 
These inequalities are only violated for ququarts (two of 
the inequalities) and eight-level quantum systems (the 
third inequality), but not for qutrits. Another example 
of a tight inequality is the noncontextuality inequality for 
qutrits of Klyachko et al. Q , which indeed was derived 
by means of the correlation polytope method. However, 
this latter inequality does not have a state-independent 
quantum violation. 

Obtaining all tight inequalities is, in general, a hard 
task. The correlation polytope is characterized by the 
number of settings and outcomes of the considered sce- 
nario. While there are algorithms that find all the facets 
of a given polytope, the time required to compute them 
grows exponentially as the number of settings increases. 
Therefore, this method can only be applied to simple 
cases with a reduced number of settings @, H, 13 1. Given 
the facets of the polytope, in a next step one can try to 
find quantum observables that exhibit a maximal gap be- 
tween the maximal noncontextual value and the quantum 
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prediction. 

In this paper we approach the problem differently. For 
many situations, the quantum observables are already 
known, and it remains to find inequalities that are tight 
and optimal and, in addition, may exhibit SIC. Thus we 
first describe the noncontextuality polytope for a given 
set of observables and a given list of admissible contexts. 
A noncontextuality inequality is then an affine hyper- 
plane that does not intersect this polytope. We then in- 
troduce a method for maximizing the state-independent 
quantum violation via linear programming. The resulting 
linear program can be solved with standard optimization 
routines, and the optimality of the solution is guaran- 
teed. As an application we derive the optimal inequality 
for several state independent scenarios, in particular an- 
alyzing a recently discovered qutrit scenario [14j . Using 
our method, we find noncontextuality inequalities with 
state- independent violation and the fewest number of ob- 
servables and contexts. These inequalities turn out to 
be in addition tight and hence provide the most funda- 
mental examples of inequalities with state-independent 
violation. 

Contextuality scenarios, the noncontextuality poly- 
tope, and noncontextuality inequalities. — We start from 
some given dichotomic [15j quantum observables 
A\, A%, . . . , A n . A context c is then a set of indices, such 
that Ak and An are compatible whenever k,£ 6 c, i.e., 
[j4fc,^] — 0. For example if A\ and A2 are compatible, 
then valid contexts would be {1}, {2}, and {1, 2}. As we 
see below, it may be interesting to consider only a cer- 
tain admissible subset £ of the set of all possible contexts 
{c}. The observables Ai, . . . ,A n , together with the list 
of admissible contexts £, form the contextuality scenario. 

The set of all (contextual as well as noncontextual) 
correlations for such a scenario can be represented by the 
following standard construction. We first use that, if Ah 
and Ag are compatible, then the expectation value of Ak 
is not changed whether or not Ag is measured in the same 
context. Thus, instead of considering all correlations, it 
suffices to only consider the vector v — (v £ | c G €), 
where v c is the expectation value of the product of the 
values of the observables indexed by c. For example, 
for the contexts {1}, {2}, {1,2}, a contextual HV model 
may with equal probability assign the values {+1}, {+!}, 
{— 1,+1}, or { — 1}, { — 1}, {+1,-1}, respectively, yield- 
ing v = { Vl ,v 2 , ui, 2 ) = (1/2, 1/2, -1). 

In the simplest noncontextual HV model, however, 
each observable has a fixed assignment a = (01, . . . , a n ) £ 
{ — 1, l} n for the observables A\, . . . ,A n , and accordingly 
each entry in v is exactly the product of the assigned val- 
ues, i.e., v c = Ylkec a k- The most general noncontextual 
HV model predicts fixed assignments ay' with proba- 
bilities pi, and hence the set of correlations that can be 
explained by a noncontextual HV models is characterized 
by the convex hull of the models with fixed assignments, 



thus forming the noncontextuality polytope. 

Then, a noncontextuality inequality is an affine bound 
on the noncontextuality polytope, i.e., a real vector A 
such that 77 > A • v for all correlation vectors v that orig- 
inate from a noncontextual model: 
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for any assignment a = (a±, . . . , a n ) £ { — 1, 1} U . 

In quantum mechanics, in contrast, the measurement 
of the entry v c corresponds to the expectation value 
(J| fc6c Afc) p , where p specifies the quantum state. Thus 

the value of X-v predicted by quantum mechanics is given 
by (T(A)) P , with 



T(A) = ^A,n^- 



(2) 
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If the expectation value exceeds the noncontextual limit 
77, then the inequality demonstrates contextual behavior, 
yielding the quantum violation 



v _ max, (T(A)) p ^ 
V 



(3) 



An inequality is optimal, if the violation is maximal for 
the given contextuality scenario. In general, this opti- 
mization is difficult to perform and it is not always clear 
that an optimal inequality also yields the most significant 
violation (l6| . 

Optimal state independent violation and tight inequal- 
ities. — However, if we require a state independent viola- 
tion of the inequality, without loss of generality, T(A) = 1 
and hence the optimization over the quantum state g 
vanishes. Then, the coefficient vector A and the noncon- 
textuality bound 77 are optimal if 77 is minimal under the 
constraint T(A) = 1 and if the inequalities in Eq. ([1]) are 
satisfied. That is, we ask for a solution (77*, A*) of the 
optimization problem 



minimize: 



subject to: T(A) = 1 and 

Eq. U) holds for all a. 



(4) 



This optimization problem is a linear program and such 
programs can be solved efficiently by standard numer- 
ical techniques and optimality is then guaranteed. We 
implemented this optimization using CVXOPT 17[ for 
Python, which allows us to study inequalities with up to 
n = 21 observables and |£| = 131 contexts. Note that 
this program also solves the feasibility problem, whether 
a contextuality scenario exhibits SIC at all. This is the 
case, if and only if the program finds a solution with 
77 < 1 and thus V > 0. 

The optimal coefficients A* are, in general, not unique 
but rather form a polytope defined by Eq. ([T]) with 77 = 
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vi = (1,0,0) u 5 = (1,0,-1 
v 2 = (0,1,0) v e = (1,-1,0) v B 
v 3 = (0,0,1) ur = (0,1,1) «o 

«4 = (0,1,-1) « 8 = (1,0,1) 
«9 = (1,1,0) 



(-1,1,1) 
(1,-1,1) 
(1,1,-1) 
»D = (1,1,1) 



FIG. 1. Graph of the compatibility relations between the 
observables for the Yu-Oh scenario. Nodes represent vectors 
\vi) [or the observables A% defined in l[5]l] and edges represent 
orthogonality (or compatibility) relations. 
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TABLE I. Coefficients X L of inequalities for the Yu-Oh sce- 
nario. The column c labels the different contexts, YO the co- 
efficients in the inequality of Ref. [TJ[ , opt 2 an optimal tight 
inequality with contexts of maximal size 2, opt 3 an optimal 
tight inequality with contexts of all sizes. For compactness, 
the coefficients in the column YO have been multiplied by 
50/3, for the column opt 2 by 52/3 and for the column opt 3 
by 83/3. The row labeled "A— 15" shows the coefficients for 
the contexts {A}, {B}, {C}, {D} and the row labeled "*,A- 
D" shows the coefficients for {4, A}, {8, A}, {9, A}, {5,B}, 
{7,B}, {9,B}, {6,C}, {7,C}, {8,C}, {4, D}, {5,D}, {6,D}. 



rf . This leaves the possibility to find optimal inequalities 
with further special properties. There are at least two im- 
portant properties that one may ask for. Firstly, from an 
experimental point of view, it would be desirable to have 
some of the coefficients X £ = 0, since then the context c 
does not need to be measured. In general, it will depend 
on the experimental setup, which coefficients A c = yield 
the greatest advantage. For the sequential measurement 
schemes it is natural to choose the longest measurement 
sequences. Secondly, there might be tight inequalities 
among the optimal solutions: An inequality is tight, if 
the affine hyperplane given by the solutions of r\ = A • x 
is tangent to a facet of the noncontextuality polytope. 
This property can be readily checked using Pitowsky's 
construction [3J : Denote by p the affine dimension of the 
noncontextuality polytope and choose those assignments 
a, for which Eq. (JTJ is saturated. Then, the inequal- 
ity is tangent to a facet if and only if the affine space 
spanned by the vertices v = (Y\ kec ak | c G £) is (p — 1)- 
dimensional. 

Furthermore, we mention that the condition of state 
independence might be loosened to only require that the 
quantum violation is at least V for all quantum states. 
This corresponds to replacing the condition T(A) = 1 
by the condition that T(A) — 1 is positive semidefinitc. 
Then, the linear program in Eq. (jj) becomes a semidef- 
inite program, which still can be solved by standard nu- 
merical methods with optimality guaranteed. However, 
for the examples that we consider in the following, the 
semidefinite and the linear program yield the same re- 
sults. 

Most fundamental noncontextuality inequalities. — We 
now apply our method to the SIC scenario for a qutrit 
system introduced by Yu and Oh [14| . Qutrit systems are 



of fundamental interest, since no smaller quantum system 
can exhibit a contextual behavior Q . It has been shown, 
that this scenario is the simplest possible SIC scenario 
for a qutrit 18 1. 

For a qutrit system, the dichotomic observables are of 
the form 



Ai = t-2\vi)( Vi \. 



(5) 



In the Yu-Oh scenario, there are 13 observables defined 
by the 13 unit vectors \vi) provided in Fig. [TJ In the 
according graph, each operator is represented by node 
i G V of the graph G = (V, E) and an edge G E 

indicates that \vi) and \vj) are orthogonal, (vj\vi) — 0, so 
that Ai and Aj are compatible. The original inequality 
takes into account all contexts of size one and two, £yo = 
{{1}, . . . , {D}} U E and the coefficients were chosen to 
A £ = -3/50 if C G E and X £ = 6/50 else. This yields an 
inequality with a state-independent quantum violation of 
V = 1/24 » 4.2%. 

With the linear program we find that the maximal vi- 
olation for the contexts Cyo is V = 1/12 ss 8.3% and 
thus twice that of the inequality in Ref Interest- 
ingly, among the optimal coefficients A* there is a solu- 
tion which is tight and for which the coefficients A4J van- 
ishes, cf. Table U column "opt2" for the list coefficients. 
We find that up to symmetries, A^ is the only context 
that can be omitted while still preserving optimality. 

In order to demonstrate the practical advantage, let 
us discuss the recent experimental values obtained for 
the Yu-Oh scenario [3, FIG. 2]. For those values, the 
original Yu-Oh inequality is violated by 3.7 standard de- 
viations. But if the same data is evaluated using our 
optimal inequality "opt2" , the violation increases to 7.5 
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standard deviations. We mention, however, that the par- 
ticular experimental setup implements the same observ- 
able in different context differently, thus easily allowing 
a noncontextual HV model explaining the data l20ll. A 
setup avoiding such problems is described in Ref. [21] . 

The maximal contexts in the Yu-Oh scenario are of 
size three, and hence it is possible to include also the 
corresponding terms in the inequality, i.e., we extend the 
contexts £yo by the contexts {1, 2, 3}, {1, 4, 7}, {2, 5, 8}, 
and {3,6,9}. Since this increases the number of param- 
eters in the inequality, there is a chance that this case 
allows an even higher violation. In fact, the maximal vi- 
olation is V = 8/75 10.7%. Again, it is possible to 
find tight inequalities with vanishing coefficients, and in 
particular the context {1, 2, 3} can be omitted; the list of 
coefficients is given in Table [H column "opts" . 

Further examples. — Our method is applicable to all 
SIC scenarios, providing the optimal inequality. We men- 
tion two further examples: (i) The "extended Peres- 
Mermin square" uses as observables all 15 products of 
Pauli operators on a two-qubit system, (a^ <g> a u ) [22j |. 
The optimal violation is V = 2/3, where only contexts 
of size three need to be measured and A £ = 1/15, ex- 



cept A. 



xx,yy,zz 



yx,zy 



= A 



xy,yz,zx 



= —1/15. Among 



the optimal solutions no simpler inequality exists, (ii) 
The 18 vector proof 23] of the Kochen-Specker theorem 
uses a ququart system and 18 observables of the form 
([5]). For contexts up to size 2 the maximal violation is 
V = 1/17 « 5.9% (cf. [H), while including all context 
the maximal violation is V = 2/7 w 28.6% (cf. 5]). The 
situation where only contexts up to size 3 are admissi- 
ble has not yet been studied and we find numerically a 
maximal violation of V « 14.3%. 

Conclusions. — Contextuality is suspected to be one 
of the fundamental phenomena in quantum mechanics. 
While it can be seen as the underlying property of the 
nonlocal behavior of quantum mechanics, so far no meth- 
ods for a systematic investigation have been developed. 
We here showed that Pitowsky's polytope naturally gen- 
eralizes to the noncontextual scenario and hence the 
question of a full characterization of this noncontextu- 
ality polytope arises. This can be done via the so-called 
tight inequalities. On the other hand, among the most 
striking aspects where contextuality is more general than 
nonlocality is that the former can be found to be inde- 
pendent of the quantum state. For this state- independent 
scenario, we showed that the search for the optimal in- 
equality reduces to a linear program, which can be solved 
numerically with optimality guaranteed. We studied sev- 
eral cases of this optimization and find that in all those 
instances one can construct noncontextuality inequalities 
with a state independent violation that are, in addition, 
tight. This is in particular the case for the most funda- 
mental scenario of state independent contextuality 14 1 



size three. We hence lifted the Yu-Oh scenario to the 
same fundamental status as the CHSH Bell inequality 

25], which is the simplest scenario for nonlocality. Our 
state-independent tight inequalities are particularly suit- 
able for experimental tests and hence we expect that they 
stimulate experiments to finally observe SIC in qutrits 

21[. 
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and we presented two essentially different inequalities — 
one involves at most contexts of size two, the other of 
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